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ABSTRACT 

Using the 5 year WMAP data, we re-investigate claims of non-Gaussianities and asymmetries de- 
tected in local curvature statistics of the 1 year WMAP data. In Hansen et al 2004, it was found that 
the northern ecliptic hemisphere was non-Gaussian at the ^ 1% level testing the densities of hill-, lake 
and saddle points based on the second derivatives of the CMB temperature map. The 5 year WMAP 
data has a much lower noise level and better control of systematics. Using these, we find that the 
anomalies are still present at a consistent level. Also the direction of maximum non-Gaussianity re- 
mains. Due to limited availability of computer resources, Hansen et al. 2004 were unable to calculate 
the full covariance matrix for the test used. Here we apply the full covariance matrix instead of 
the diagonal approximation and find that the non-Gaussianities disappear and there is no preferred 
non-Gaussian direction. We compare with simulations of weak lensing to see if this may cause the 
observed non-Gaussianity when using diagonal covariance matrix. We conclude that weak lensing 
does not produce non-Gaussianity in the local curvature statistics at the scales investigated in this 
paper. The cause of the non-Gaussian detection in the case of a diagonal matrix remains unclear. 
Subject headings: cosmic microwave background — cosmology: observations — methods: statistical 



1. INTRODUCTION 

During recent years, the investigations of the cos- 
mic microwave background (CMB) has proved to 
be the most compelling addition to our understand- 
ing of the early universe. Observations of the 
CMB anisotropies, like those obtained by the Wilkin- 
son Microwave Aniso t ropy Probe (WMAP ) experiment 
([Bennett et all 120031: iHinshaw etall I2007D . have pro- 
vided us with profound insight on the composition of 
structure in our universe. Combined with previous ex- 
perimental knowledge and a sound theoretical frame- 
work, the concordance model of ACDM has been estab- 
lished. 

The ACDM model relies on the framework of infia- 
tion. Infiation was initially pro posed as a solutio n to the 
horizon and flatness problem (jGuth et al 1119811 ). Addi- 
tionally, it established a highly successful scenario for the 
formation of primordial density perturbations, providing 
the required seeds for the large-scale structures in the 
universe. Eventually, these later gave rise to the temper- 
ature anisotropies in the cosmic microwave background 
radiation that we observe today. It is assumed that Gaus- 
sian fluctuations in the vacuum field have given rise to 
these perturbations, so the fluctuations in the CMB map 
we observe today are also expected to be near Gaussian. 

Inflation predicts that the observed universe should 
be statistically isotropic on large scales. However, 
during recent years, various anomalies that contradict 
statis tical isotropy have been discovered in the WMAP 
data (Ide Oliveir a-Cost a et al.ll2004[: lEriksen et al.ll2004al: 
[Hansen et aiTBoOQ: Hoftuft et al.[ 120091: lEriksen et al.l 
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Groeneboora fc Erikseiil l2009HGroeneboom etall 12009 : 
Vielva et al.ll2004h . If these effects are confirmed by the 
data from the Planck experiment, various anisotropic 
universe models should be seriously considered. 

In order to investigate the properties of different infla- 
tionary models, it is important not only to pursue devi- 
ations from statistical isotropy, but also any deviations 
from Gaussianity. Single-field inflation models usually 
predict a CMB statistically close to Gaussian, but more 
exotic models may give rise to a larger contribution of 
non-Gaussianity. In this paper, we focus on a method 
for testing general deviations from Gaussianity that are 
not obviously connected to infiationary non-Gaussianity. 
However, some effort has been made to use the local cur- 
vature to estimate th e non-Gaussianity parameter /jvl 
(ICabella et al 1120051) . 

IDore et al.l ( 20031 ) presented a framework for investi- 
gating non-Gaussianities based on the properties of the 
local curvature of the CMB temperature map. By calcu- 
lating the second-order derivatives, it is possible to clas- 
sify each pixel as a "hill", "lake" or "saddle" based on 
the eigenvalues of the Hessian matrix. If the map is first 
smoothed with a Gaussian beam, it is possible to ex- 
tract the local curvature properties on scales given by 
the FWHM of the beam. A temperature threshold Tt 
is introduced, where CMB temperature values below a 
provided threshold are ignored. Starting with a negative 
temperature threshold, the fraction of hill, lake and sad- 
dle points that have temperatures above the threshold 
are counted. By performing this analysis on simulated 
isotropic Gaussian maps, it is possible to estimate what 
the fraction of hills, lakes and saddle points should be for 
increasing Tt. This graph is then compared with a similar 
examination of experimental data using a test. It is 
is also possible to pr edict this plot theoretically, as done 
bv IDore et al.l (|2003D . However, several issues will com- 
plicate the task. First, one needs to apply a galaxy mask 
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that effectively removes about 20% of the data. There 
are also point sources to be considered, along with the 
possibility of excluding large chunks of the sky in order 
t o focus the ana l ysis o n specific areas. 

[Hansen et al.l ( 20041 ) applied their framework on the 
1-yr WMAP data, where no statistical deviation from 
Gaussianity was found on a full-sky coverage. However, 
when the authors analyzed independent hemispheres on 
the CMB sky, they found a non-Gaussian signature on 
hemispheres centered near the ecliptic poles. The au- 
thors therefore suggested that the effect might be due 
to systematics, an d that the direction s eems consistent 
with the results bv lEriksen et al.l ()2004aD . However, due 
to limited computer resources, the authors built their 
statistics using only 128 and 512 Gaussian simulations, 
which may be too few for obtaining convergence. For the 
same reason, they were unable to obtain a full covcrgcd 
covariance matrix for the analysis and they therefore 
ignored the correlations between threshold levels. In ad- 
dition, the 1-year WMAP data are much more noisy com- 
pared to the 5-year data. In this paper, we re-analyze the 
5-year WMAP data using an independent and more ro- 
bust code with a greatly increased number of Gaussian 
simulations. We also test whether weak lensing may con- 
tribute to a non-Gaussian signal in the local curvature 
statistics. 

2. METHOD 

The methods in this p aper are based on the w ork done 
bv [Hansen et al.l (|2004f) and iDore et"all (|2003[ ). We are 
interested in studying the curvature properties of the 
temperature fluctuations in the CMB map. Every point 
on a two-dimensional surface embedded in 3-space can 
be characterized as being a hill, lake or saddle. In the 
case of the CMB, we usually express coordinates on §^ 
in spherical coordinates (6*, i^). A hill/lake is where the 
curvature is negative/positive in both the 6 and (p direc- 
tions, whereas a saddle point has negative curvature for 
one coordinate and positive for the other. While increas- 
ing the threshold value Tt, we count the fraction of hill, 
lake and saddle points for all the pixels with temperature 
larger than Tt- A simulated CMB map with the classified 
curvature areas added on top for illustrative purposes is 
depicted in Figure [1] 

2.1. The local curvature 

The CMB map is first smoothed with a beam corre- 
sponding to the scales that we are interested in. We then 
normalize the temperature map T with its standard devi- 
ation and calculate its various first and second-order co- 
variant derivatives given the spherical coordinates 6 and 
</.. That is, we obtain ^, f , ^, and For 
each pixel in the the normalized temperature map T, we 
calculate the Hessian matrix and calculate its eigenva l- 
ues Ai and A2 (as described in lMonteserin et all ()2005() ). 
Then, depending on the sign of the eigenvalues, we clas- 
sify every point on the sphere as a hill (Ai > 0, A2 > 0), 
lake (Ai < 0, A2 < 0) or saddle point (Ai,A2 < 0). 

By performing this analysis on simulated isotropic 
Gaussian maps, it is possible to estimate what the mean 
fraction of hills, lakes and saddle points should be for 
increasing Tt, together with the standard deviation. An 
example of such a graph is presented in Figure [2l 




Fig. 1. — The classification of liills (red parts), lakes (blue parts) 
and saddles (neutral, in between). The original CMB map is added 
in the background for illustrative purposes. The beam FWHM used 
in this example is 500' on an isotropic realization of the best-fit 
ACDM power spectrum. 




Temperature threshold 

Fig. 2. — The fraction of hills, lakes and saddle pixels in a simu- 
lated Gaussian map with nai^c = 512. The simulation is generated 
from the best-fit powerspectrum from the WMAP 5 year release. 
The black line depicts the mean, while the shaded area is the 68% 
confidence area. 

2.2. The extended mask 

When working with experimental CMB data, we need 
to remove a portion of the sky due to foreground residu- 
als. In this paper, we operate with s everal ma sks, most 
notably the WMAP KQ85 sky cut ()Gold et a l. 2008), 
which removes 18% of the sky. The edges of the mask 
need special consideration, due to the differentiation of 
the map. We therefore perform the following procedure 
for generating suitable CMB masks: 

1. Start with the 5-year WMAP KQ85 or KQ75 mask 
without point source holes. 

2. The mask is expanded by b°, where b is the 
smoothing angle FWHM to avoid effects of the 
s moothed mask. Note that this step was not used 
in [Hansen et al.l (l2004f ). We thus use an extended 
mask which is considerably much larger than the 
one used in [Hansen et al., (,20041 ). As pointed out 



3 




Fig. 3. — The KQ85 mask expanded with 5° and point-sources 
added. The expansion from differentiation is indicated in grey. 

below, we have compared results also with the 
smaller extended mask for which this step was 
omitted. 

3. Point source holes are added to the mask after ex- 
pansion. If added before smoothing, the smeared 
point sources would nearly fill the mask. 

4. The mask is differentiated, and then normalized. 
By normalized, we mean that each pixel pi is set 

to = Jf^. 

5. If a pixel in a normalized, differentiated map has 
value above a certain threshold, the pixel is masked 
out. We have used the threshold 0.02 throughout 
this paper. 

2.3. statistics 

We now present a method for estimating how much a 
CMB data set deviates from Gaussianity, based on stan- 
dard statistical tests. We calculate the covariance matrix 
C from the fraction of hills, lakes and saddles for 40000 
Gaussian simulations, 
on the form 



We implement a standard x -test 



= d^c-M, (1) 

where d = — {X)Tt and is a hill, lake or sad- 
dle density at threshold Tt- We then estimate the on 
10000 simulations, and obtain a histogram of the distri- 
bution, as presented in Figure |4l When performing the 
analysis on WMAP data, we first calculate the x^ and 
then compare it to the pre-calculated x^ distribution. 
We then count the percentage of the x^s that are above 
the experimental, i.e. if 32% of the simulated x^ values 
are larger than the x^ from experimental data this cor- 
responds to 1(7 deviation from the Gaussian expectation. 
A value of 5% would be consistent with a 2a deviation 
from the Gaussianity. 

As an alternative method we also use a diagonal covari- 
ance matrix for calculating the values, thus ignoring 
correlatio ns between threshold levels. This is the method 
applied in [Hansen et al.l (|2004[ ). 

3. DATA 

We consider the publicly available 5-year WMAP data 
(jBennett et al.l l2003t iHinshaw etall I2007D that can be 
obtained from the LAMBDA0 site. We primarily per- 
form the analyses on a combined V and W frequency 

* http://lambda.gsfc.nasa.gov/ 
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X from data 



Fig. 4. — An example of a table with obtained from Gaus- 
sian simulations (histogram) and a from real data (line). The 
fraction of simulated x^ above the x^ from real data describes its 
deviation from the Gaussian mean. 

band (94 and 61 GHz) in order to keep noise and fore- 
grounds to a minimum. In addition to the 5-yr WMAP 
data, we also simulate 50000 Gaussian CMB maps based 
on the best-fit angular powerspectrum from the 5-year 
WMAP release. We convolve the simulated maps with 
the combined V and W instrumental beam before adding 
Gaussian noise. All CMB maps considered in this paper 
have Healpix resolution iVsidc = 512. 

We focus on several different divisions of the sky. First, 
wc consider the full-sky WMAP data including the ex- 
tended KQ85 and the KQ75 mask from Section [2?2] We 
then investigate the northern and southern galactic hemi- 
spheres individually, using only the KQ85 mask. Due 
to the alignment of the WMAP satellite with the eclip- 
tic plane, we also include an analysis of the northern 
and southern ecliptic hemispheres. We then perform the 
full-sky analysis on a co-added Q-band data set, a co- 
added V-band set and a co-added W-band set. Finally, 
we re-analyze the co-added V-l-W data sets on hemi- 
spheres centered around all the pixel centers on a map 
with iVsido = 2. From these results, we obtain maps with 
preferred directions depicting the amount of deviation 
from Gaussianity on the sky. 

4. RESULTS 

Here we present the results, first with a full covariance 
matrix and then without. 

4.1. Results with full covariance matrix 

The results from the V-l-W full sky and north/south 
galactic/ecliptic analyses are presented in Table [T] The 
most interesting result from Table [T] is the ~ 1% devia- 
tion from Gaussianity in the full-sky at 3° when counting 
saddles. A general trait throughout the analysis is that 
scales above and below 3° have deviation from Gaussian- 
ity no more than 2a. We continue by analyzing whether 
the 1% detection on the full-sky can be observed in the 
various WMAP frequency bands alone. We therefore 
perform the same full-sky analysis on the Q, V and W 
frequency bands. The results are shown in Table [3l and 
are consistent with the combined V+W full-sky analy- 
sis, where a 2a deviation from Gaussianity is visible in 
all band for saddles at 3°. 

Mot ivated by the asymmetry found in [Hansen et al.l 
(|2004D . we perform the analysis on all hemispheres cen- 
tered around pixels on a map with iVside = 2 using the 
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Fig. 5. — The results from an non-Gaussian analysis of hemispheres centered around pixels on a map with A'sidc = 2. Darker pixels 
correspond to higher deviations from Gaussianity on a hemisphere centered around this pixel. The values of the pixels indicate the 
percentage of simulations with a higher x'^- Whiter pixels correspond to < la levels. 



co-added V+W data set. This analysis used the ex- 
tended KQ85 mask. The results are shown in Figure 
[5l where darker pixels represent larger deviations from 
Gaussianity on a hemisphere centered around the pixel. 
Here, whiter pixels correspond to < la levels. These re- 
sults agree with the full-sky analysis and with the eclip- 
tic/galactic hemispherical analyses, and show that there 
is no apparent preferred direction in either of the data 
sets. 

Finally, we perform a test of the combined results. 
We construct a data vector containing all fractions of 
hills and lakes for all thresholds and smoothing scales 
(1°, 3° and 5°). The covariance matrix including all cor- 
relations between hills and lakes and between different 
scales, is calculated from Gaussian simulations. We esti- 
mate a value for both WMAP data and simulations. 
The final comparison shows that 80% of the simulations 
have a higher value than the simulations for the com- 
bined V -t- W channel. In other words, when using the 
full covariance matrix, there is full consistency with the 
Gaussian hypothesis. 

4.2. Comparison with earlier work 

[Hansen et al.l ()2004l) found very little evidence for de- 
viation from Gaussianity in the 1-year WMAP data 
when analyzing full-sky CMB maps using the Kp2 galaxy 
mask. This is in good agreement with the findings in 
this paper, where wc have shown there is at most a 2a 
deviation from Gaussianity in the saddles at 3° in the 
combined V-l-W data. We stress that while the num- 
ber o£_simulated Gaussian maps used for the test 
by [Hansen ct aP (2004') was 512, we have used 50 000. 
Where .Hansen et al., (,2004,1 employed a covariance ma- 



trix that was strictly diagonal, we have operated with a 
full covariance matrix including correlations between dif- 
ferent thresholds. In addition, the 1-year WMAP data 
are in general more contaminated by instrumental noise 
than the 5-year data. 

When using the full covariance matrix, the directional 
analysis show that there isn't any preferred direction for 
deviations from Gaussianity. Figure [S] shows that the 
directions for various scales seem to be randomly scat- 
tered, and the significance is very low. In fact, only at 
3° there seems to be a 2a hint, but any exact directions 
are n on-existent. This is i n disagreement with the results 
from [Hansen et all (|2004[ ). who claimed that there exists 
a maximum of no n- Gaussianity on hemispheres centered 
at the ecliptic poles. 

However, when we ignore correlations and use a strictly 
diagonal covariance matri x, the results are more in agree- 
ment with Hansen et al] ()2004[ ). These results are pre- 
sented in Table |2l where we note several deviations at 
the less than 1% level. Also, when performing the direc- 
tional analysis, we find that the non-Gaussian signal has 
a clear maximum close to (but not directly on) the eclip- 
tic north-pole, as seen in figure [5] and Tabl ed This is in 
agreem ent with the direction described by [Hansen et al.1 
(|2004f ) . As n oted above, here wc use a larger extended 
mask than in [Hansen et al.l (|2004[ ). We have also tested 
with the smaller extended mask and find results simiar 
to the ones presented in Table [H We conclude that by 
using a diagonal approximation to the correlation ma- 
trix, we still find non-Gaussianities and asymmetries over 
several scales and with different masks, including the full 
covariance matrix lowers the effect of the non-Gaussian 
signatures and the hemispherical anisotropy. 
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Fig. 6. — The results from a non-Gaussian analysis using a diagonal covariance matrix, of hemispheres centered around pixels on a map 
with A^sidc = 2. Darker pixels correspond to higher deviations from Gaussianity on a hemisphere centered around this pixel. The values of 
the pixels indicate the percentage of simulations with a higher x^- Whiter pixels correspond to < Icr levels. 



4.3. Tests on weak tensing simulations 

In order to investigate whether the deviations may be 
a resuh of weak Icnsing in the CMB, wc performed an- 
other te s t. W e used the freely available LensPijfl code 
([Lewis I (120051 )) to simulate 300 Gaussian CMB maps, 
and 300 CMB maps with weak Icnsing. For the simula- 
tions we used the best-fit power spectrum and Icnsing po- 
tential provided by the WMAP team. The procedure of 
counting hills, lakes and saddles as described above was 
applied to all the simulations. From the Gaussian sim- 
ulations we found mean values and standard deviations 
of hill, lake and saddle fractions at different smoothing 
scales. We then compared these with the results from 
the simulations with weak lensing. We found no sig- 
nificant deviation in the local curvature of the simula- 
tions with lensing. We therfore conclude that the devi- 
ations from n o n- Gau ssianity described in this paper and 
Hansen et al.l (|2004[ ) are not caused by weak lensing in 
the CMB. 

5. CONCLUSIONS 

We have developed an independent framework for esti- 
mating deviations from Gaussia nity in CMB data based 
on the methods est ablished by iDore et all ()2003[ ) and 
Hansen et al.l (|2004f ). The methods used are model- 
independent, and do not share any obvious connections 
with non-Gaussianity frameworks of known physical ori- 
gin. By counting the fraction of lakes, hills and saddles 
in simulated Gaussian maps while increasing the tem- 
perature threshold, we have built a distribution for what 
is expected for Gaussian maps. We then compared ex- 

^ http:/ /cosmologist.info/lcnspix/ 



perimental data to this distribution, determining the de- 
viation from the Gaussian assumption. Wc then con- 
sidered a combined V -I- W full-sky data set with the 
extended KQ85 and KQ75 mask, and found evidence of 
a ~ 1% deviation from Gaussianity on scales around 3°. 
We also analyzed other scales as well as the north and 
south galactic/ecliptic hemispheres seperately, but dis- 
covered no deviation from Gaussianity greater than 2a. 
We continued by performing an analysis on each of the 
hemispheres centered around a pixel on a Healpix map 
with iVsidc = 2 using the combined V+W data. We pro- 
duced directional maps for hills, lakes, saddles using the 
three scales, and found no evidence for a preferred di- 
rection in either of the maps. Finally, we calculated the 
combined from all our results, which resulted in an 
overall agreement with Gaussianity. We conclude that 
there is no significant evidence for non-Gaussianities or 
asymmetries i n the WMAP data based on this test. 

However, in Hansen et al.l (j2004[) . it was found that the 
northern ecliptic hemisphere was non-Gaussian based on 
similar local curvature measurements. There was how- 
ever one large difference in the method used: In that 
work, a diagonal approximation to the covariance matrix 
was applied. Repeating our analysis with a diagonal co- 
varian ce matrix we obtain similar results as Hansen et al.l 
(j2004D . Taking into account correlations between thresh- 
olds, the non-Gaussianity disappears. We have compared 
Gaussian CMB simulations to CMB simulations with 
weak Icnsing to sec whether the hill, lake or saddle densi- 
ties arc different in the two sets of simulations. No signif- 
icant differences were found. We therefore conclude that 
weak lensing may not cause the non-gaussianity found 
when using a diagonal correlation matrix. It is still un- 
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TABLE 1 

Summary of deviations from Gaussianity in the combined 
V+W data 



Area/Scale 


60' 


180' 


300' 






Full sky H 


h KQ85 mask 




hills 


27% 


6.7% 


3.87% 


saddles 


69% 


1.13% 


11% 


lakes 


46% 


15% 


6.5% 




Full sky H 


h KQ75 mask 




hills 


48% 


10% 


14% 


saddles 


85% 


1.33% 


7% 


lakes 


92% 


14% 


16% 


Northern galactic hemisphere 


hills 


86% 


30% 


46% 


saddles 


75% 


13% 


25% 


lakes 


73% 


18% 


38% 


Southern galactic hemisphere 


hills 


39% 


5% 


10% 


saddles 


78% 


5% 


54% 


lakes 


53% 


5% 


12% 


Northern ecliptic hemisphere 


hills 


64% 


28% 


15% 


saddles 


83% 


5% 


20% 


lakes 


69% 


45% 


25% 


Southern ecliptic hemisphere 


hills 


73% 


13% 


13% 


saddles 


75% 


15% 


54% 


lakes 


62% 


14% 


17% 



Note. — The values correspond to the fraction of x'^ from 
isotropic simulations that have a higher than the analysis of the 
data set. A low percentage would thus indicate a non-Gaussianity. 



TABLE 2 

Summary' of deviations from Gaussianity in the combined 
v+w data usinc; a strictly diagonal cmvariance matrix 



Area/Scale 


60' 


180' 


300' 




Full sky + KQ85 mask 


hills 


12% 


0.95% 


0.68% 


saddles 


11% 


0.034% 


10% 


lakes 


11% 


0.67% 


0.64% 




Northern ecliptic 


hemisphere 




hills 


13% 


0.38% 


4.2% 


saddles 


11% 


0.34% 


11% 


lakes 


3.8% 


0.67% 


4.2% 




Southern ecliptic 


hemisphere 




hills 


39% 


30% 


15% 


saddles 


78% 


4.0% 


87% 


lakes 


29% 


3.4% 


12% 



Note. — The values correspond to the fraction of from 
isotropic simulations that have a higher than the analysis of the 
data set. A low percentage would thus indicate a non-Gaussianity. 



TABLE 3 

Summary of deviations from Gaussianity in full-sky 
combined bands 



Area/Scale 


60' 180' 


300' 




Full-sky combined W-band 


hills 

saddles 

lakes 


20% 7.5% 
66% 2.8% 

37% 27% 


5.9% 

24% 
5% 


Full-sky combined V-band 


hills 

saddles 

lakes 


24% 11% 
39% 1.61% 
65% 27% 


8.2% 
13% 
7% 


Full-sky combined Q-band 


hills 

saddles 

lakes 


15% 6.6% 
76% 1.73% 

34% 25% 


3.13% 
8.9% 
5% 



Note. — The values correspond to the fraction of from 
isotropic simulations that have a higher x^ than the analysis of the 
data set. A low percentage would thus indicate a non-Gaussianity. 



clear what causes the detection of non-Gaussianity when 
a diagonal correlation matrix is used. Even though the 
test is not optimal when correlations are ignored, we 
are still comparing the data to simulations for which an 
identical procedure (i.e. diagonal approximation to the 
covariance matrix) has been applied. If the best fit model 
estimated from the data is correct, one would expect sim- 
ulations based on this model to have the same statistical 
properties as the data, no matter which statistical test is 
performed. There is thus still a discrepancy between data 
and simulations based on the model which best fits the 
data. Whether this discrepancy is a statistical fluke or 
may arise from systematic errors/cosmology and whether 
it is related to other asymmetries is still unclear. 
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